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Abstract. Accuracy levels of metres per second require the fundamental concept of 'radial velocity' for stars 
and other distant objects to be examined, both as a physical velocity, and as measured by spectroscopic and 
astrometric techniques. Already in a classical (non-relativistic) framework the line-of-sight velocity component 
is an ambiguous concept, depending on whether, e.g., the time of light emission (at the object) or that of light 
detection (by the observer) is used for recording the time coordinate. Relativistic velocity effects and spectroscopic 
measurements made inside gravitational fields add further complications, causing wavelength shifts to depend, e.g., 
on the transverse velocity of the object and the gravitational potential at the source. Aiming at definitions that are 
unambiguous at accuracy levels of 1 m s -1 , we analyse different concepts of radial velocity and their interrelations. 
At this accuracy level, a strict separation must be made between the purely geometric concepts on one hand, and 
the spectroscopic measurement on the other. Among the geometric concepts we define kinematic radial velocity, 
which corresponds most closely to the 'textbook definition' of radial velocity as the line-of-sight component of 
space velocity; and astrometric radial velocity, which can be derived from astrometric observations. Consistent with 
these definitions, we propose strict definitions also of the complementary kinematic and astrometric quantities, 
namely transverse velocity and proper motion. The kinematic and astrometric radial velocities depend on the 
chosen spacetime metric, and are accurately related by simple coordinate transformations. On the other hand, 
the observational quantity that should result from accurate spectroscopic measurements is the barycentmc radial- 
velocity measure. This is independent of the metric, and to first order equals the line-of-sight velocity. However, it is 
not a physical velocity, and cannot be accurately transformed to a kinematic or astrometric radial velocity without 
additional assumptions and data in modelling the process of light emission from the source, the transmission of 
the signal through space, and its recording by the observer. For historic and practical reasons, the spectroscopic 
radial- velocity measure is expressed in velocity units as czb, where c is the speed of light and zb is the observed 
relative wavelength shift reduced to the solar-system barycentre, at an epoch equal to the barycentric time of 
light arrival. The barycentric radial-velocity measure and the astrometric radial velocity are defined by recent 
resolutions adopted by the International Astronomical Union (IAU), the motives and consequences of which are 
explained in this paper. 

Key words, techniques: radial velocities - techniques: spectroscopic - astrometry - reference systems - stars: 
kinematics - methods: data analysis 

trometric) concept of radial velocity - were adopted at 
the General Assembly of the International Astronomical 
Union (IAU), held in 2000. The purpose of this paper is 
to explain their background, the need for such definitions, 
and to elaborate on their consequences for future work. 
Thus, the paper is not about the detailed interpretation 
of observed spectral-line displacements in terms of radial 
motion, nor about the actual techniques for making such 
measurements; rather, it is the definition of 'radial veloc- 
ity' itself, as a geometric and spectroscopic concept, that 
is discussed. 

The need for a strict definition has become urgent in 
recent years as a consequence of important developments 
in the techniques for measuring stellar radial velocities, 
as well as the improved understanding of the many effects 



1. The need for stringent definitions 

Radial velocity is an omnipresent concept in astronomy, 
and a quantity whose precision of determination has im- 
proved significantly in recent years. Its meaning is gen- 
erally understood as the object's motion along the line 
of sight, a quantity normally deduced from observed 
spectral-line displacements, interpreted as Doppler shifts. 
However, despite its ubiquity, there has not existed any 
physically stringent definition of 'radial velocity' with an 
accuracy to match currently attainable measuring preci- 
sions. Two first such definitions - one for the result of 
spectroscopic observations, and one for the geometric (as- 
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that complicate their interpretation. We note in particular 
the following circumstances: 

Precision and accuracy of spectroscopic measurements: 
Spectroscopic measurement precisions are now reaching 
(and surpassing) levels of metres per second. In some ap- 
plications, such as the search for (short-period) extraso- 
lar planets or stellar oscillations, it may be sufficient to 
obtain differential measurements of wavelength shifts, in 
which case internal precision is adequate and there is no 
need for an accurate definition of the zero point. However, 
other applications might require the combination of data 
from different observatories, recorded over extended peri- 
ods of time, and thus the use of a common reference point. 
Examples could be the study of long-term variations due 
to stellar activity cycles and searches for long-period stel- 
lar companions. Such applications call for data that are 
not only precise, but also accurate, i.e., referring to some 
'absolute' scale of measurements. However, the transfer of 
high-precision measurements to absolute values was previ- 
ously not possible, partly because there has been no agree- 
ment on how to make such a transfer, or even on which 
physical quantity to transfer. 

In the past, a classical accuracy achieved for measur- 
ing stellar radial velocities has been perhaps 1 km s _1 , 
at which level most of these issues did not arise, or could 
be solved by the simple use of 'standard stars'. With cur- 
rent methods and instrumentation, the accuracy by which 
measured stellar wavelengths can be related to absolute 
numbers is largely set by the laboratory sources used for 
spectrometer calibration (lines from iodine cells, lasers, 
etc.). An accuracy level of about 10 m s _1 now seems 
reachable. Since any fundamental definition should be at 
least some order of magnitude better than current per- 
formance, the accuracy goal for the definition was set to 
I m s _1 . This necessitates a stringent treatment of the 
radial- velocity concept. 

Ambiguity of classical concepts: A closer inspection 
even of the classical (non-relativistic) concepts of radial 
velocity reveals that these are ambiguous at second order 
in velocity relative to the speed of light. For instance, if 
radial velocity is defined as the rate of change in distance, 
one may ask whether the derivative should be with respect 
to the time of light emission at the object, or of light re- 
ception at the observer. Depending on such conventions, 
differences exceeding I m s _1 would be found already for 
normal stellar velocities. 

Intrinsic stellar spectroscopic effects: On accuracy lev- 
els below ~ 1 km s _1 , spectral lines in stars and other ob- 
jects are generally asymmetric and shifted in wavelength 
relative to the positions expected from a Doppler shift 
caused by the motion of their centres-of-mass. Such ef- 
fects are caused e.g. by convective motions in the stel- 
lar atmosphere, gravitational redshift, pressure shifts, and 
asymmetric emission and/or absorption components. As a 
consequence, the measured wavelengths do not correspond 
to the precise centre-of-mass motion of the star. 

Relation between Doppler shift and velocity: Even if 
we agree to express the observed wavelength shift (what- 



ever its origin) as a velocity, it is not obvious how that 
transformation should be made. Should it use the classi- 
cal formula v r = cz (where c is the speed of light and 
z = (A bs — Ai a b)/Ai a b the dimensionless spectral shift), 
or the relativistic version (in which case the transverse 
velocity must either be known or assumed to be neg- 
ligible)? Differences are of second order in z, thus ex- 
ceeding 1 m s _1 already for 'normal' stellar velocities 
(^ 20 km s _1 ), and 100 m s _1 for more extreme galactic 
velocities (^ 200 km s _1 ). 

The role of standard stars: Practical radial-velocity 
measurements have traditionally relied on the use of 
standard stars to define the zero point of the velocity 
scale. While these have aimed at accuracies of the or- 
der 100 m s _1 , it has in reality been difficult to achieve 
consistency even at this level due to poorly understood 
systematic differences depending on spectral type, stellar 
rotation, instrumental resolution, correlation masks used, 
and so on. Standard stars will probably continue to play a 
role as a practical way of eliminating, to first order, such 
differences in radial-velocity surveys aiming at moderate 
accuracy. However, their relation to high-accuracy mea- 
surements needs to be clearly defined. 

Gravitational redshifts: The gravitational potential at 
the stellar surface causes all escaping photons to be red- 
shifted by an amount that varies from ~ 30 m s _1 for su- 
pergiants, ~ 30 km s _1 for white dwarfs, and much greater 
values for neutron stars and other compact objects. Even 
for a given star, the precise shift varies depending on the 
height at which the spectral lines are formed. The observed 
shift also depends on the gravitational potential at the ob- 
server, and therefore on the observer's distance from the 
Sun. 

Astrometric determination of radial motion: Current 
and expected advances in astrometry enable the accu- 
rate determin ation of stellar radial motions without using 
spectroscopy l|Dravins et alJlTQQflft . i.e., based on purely 
geometric measurements such as the secular change in 
trigonometric parallax. Comparison of such velocities with 
spectroscopic measurements could obviously give a handle 
on the intrinsic stellar effects mentioned above, but how 
should such a comparison be made? How does the astro- 
metric radial velocity differ conceptually from the spec- 
troscopically determined velocity? 

Accurate reference systems for celestial mechanics and 
astrometry: The rapid development of observational ac- 
curacies in astrometry and related disciplines has made 
it necessary to introduce new conventions and refer- 
ence systems, consiste nt with general relati vity at sub- 
microarcsecond levels l| Johnston et alJ 120001. Radial ve- 
locity, regarded as a component of space velocity, obvi- 
ously needs to be considered within the same framework. 

Cosmology: Ultimately, spectroscopic measurements of 
distant stars are also affected by cosmological redshift. 
To what extent does also the local space to nearby stars 
take part in the general expansion of the Universe? What 
is an actual 'velocity', and what is a change of spatial 
coordinates? Since such factors are generally not known 
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Fig. 1. Whether a radial velocity is measured spectroscop- 
ically or geometrically, the relevant astronomical event 
consists of several parts: the motion of the observed ob- 
ject; the emission of an electromagnetic signal from the 
object, symbolised by the open circle; the propagation of 
the signal through space; the motion of the observer; and 
the reception of the signal by the observer, symbolised 
by the black dot (figure partially adapted from iKlioneJ 
2000*3). 



of the event as illustrated in the figure. This, in turn, leads 
to the introduction of additional geometric quantities, viz. 
kinematic radial velocity and astrometric radial velocity 
(Sect. EJ. 

The remainder of this paper is organised as follows. 
Section|21contains a preliminary heuristic discussion of the 
radial- velocity concept as such; the purpose is to point out 
the inadequacy and ambiguity of standard notions, with- 
out providing a solution. Section then gives an overview 
of spectroscopic radial-velocity determinations, highlight- 
ing effects other than trivial stellar motion that influence 
the outcome of such measurements. The accurate meaning 
of the geometric and spectroscopic quantities is evaluated 
in the more technical Sections 0] and [S] leading up to the 
IAU resolutions, whose practical implications are consid- 
ered in Sect. EJ Some unsolved issues, beyond the scope 
of the present definitions, are briefly discussed in Sect.0 
The Appendix contains the full text of the two IAU reso- 
lutions. 

2. What is meant by 'radial velocity'? 

In this Section we examine the common notions of 'ra- 
dial velocity' from a somewhat naive viewpoint, in order 
to highlight some of the difficulties associated with this 
apparently simple concept. Additional complications in 
the interpretation of spectroscopic shifts are discussed in 
Sect. El 



to the spectroscopic observer, it is impossible to convert 
the observed shift into a precise kinematic quantity. 

From these examples it should be clear that the naive 
concept of radial velocity, as the line-of-sight component 
of the stellar velocity vector measured by the Doppler shift 
of the spectral lines, is far too simplistic when aiming at 
accuracies much better than 1 km s^ 1 . To arrive at a con- 
sistent set of definitions applicable to the various classes 
of observations, it is necessary to consider all the phases 
of an astronomical event (Fig. [TJ. These include the mo- 
tion of the star; the emission of a light signal from the 
star and its propagation to the observer through varying 
gravitational fields and possibly expanding space; the mo- 
tion of the observer; and the reception and measurement 
of the signal by the observer. 1 The observer has detailed 
knowledge only of the last two phases of the event (sym- 
bolised by the black dot in Fig.^J, while the interpretation 
of the previous phases requires additional assumptions or 
modelling. The result of a measurement should ideally be 
specified in a way that is neutral with respect to such in- 
terpretation. As we shall see (Sect. |SJ), this leads to the 
definition of barycentric radial-velocity measure as the de- 
sired result of a spectroscopic measurement. However, to 
relate this quantity to physical motion in even simple sit- 
uations, a model is required which incorporates all phases 



2.1. Geometric concepts 

The Encyclopedia of Astronomy and Astrophysics defines 
the radial velocity of a star as 'the com ponent of its m otion 
along the line of sight of the observer' l|Lathaml200l|) . The 
Explanatory Supplement to the Astronomical Almanac 
gives an alternative d efinition, 'the rate of change of the 
distance to an object' (|Seidelmannl ll992'). Both agree with 
common notions about radial velocity, but are they equiv- 
alent? Let us start by examining this question in a purely 
classical framework, without the complications of relativ- 
ity, but taking into account the finite speed of light (c). 

In a Euclidean metric with origin at the solar-system 
barycentre and with t denoting coordinate time, let r*(t) 
be the motion of the star, = dr*/dt its barycentric 
space velocity, r* = |r*| the barycentric distance, and u = 
r^jr* the barycentric direction to the star. Following the 
first of the two definitions quoted above, the radial velocity 
v r is the component of along u, or 

v r = u'v* , (1) 

where the prime (') denotes scalar multiplication of vec- 
tors. From r* = «r* we can write the space velocity 



du 



dr* 



u ■ 



(2) 



For conciseness, we will from here on use the word 'star' to 
denote any observed object far outside the solar system, and 
'light' to denote electromagnetic radiation from that object. u v * 



_ dr„ 

V * = ~dt dt " ' ~ dt 
Taking the scalar product with u we have, since u'u = 1 
and u'(du/dt) = ±d(u'u)/dt = 0, 

, dr* 



dt 



(3) 
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The right member apparently corresponds to the second 
definition quoted above. Comparing Eqs. and @ it 
therefore appears that we have proved the equivalence of 
the two definitions. 

However, the situation is more complex when the fi- 
nite speed of light is considered. The observation involves 
(at least) two different times, viz. the time of light emis- 
sion at the star (t*) and the time of light reception at 
the observer (t ba), cf. Fig.Q The second definition, 'the 
rate of change of the distance', is in fact ambiguous be- 
cause it is not specified which time is used to compute 
the time derivative. Clearly the t in Eq. (JSJ must be the 
same as used in describing the motion of the star, r*(t), 
which should be independent of the observer and there- 
fore corresponding to the time of light emission. However, 
when describing an observed phenomenon, such as a mea- 
surement of the line-of-sight velocity of a star, it is more 
natural to refer it to the time of light reception i b s - 2 

The two instants and i bs are related by the light- 
time equation, which for an observer at the solar-system 
barycentre (and ignoring gravitational time delay) is sim- 
ply 



(4) 



Depending on which t is used to compute the 'rate of 
change' in the second definition, we obtain by means of 
the light-time equation two different expressions for v r : 

d^obs 



dr» 
di* 

dr* 

obs 



1 



dt, 



= c 1 



dinhf 



(5) 



The difference, v' r — v" = v' r v" jc ~ v 2 /c, exceeds 1ms 1 
for \v r \ k, 20 km s" 1 and 100 m s" 1 for \v r \ k, 200 km s" 1 . 
Since relative velocities in our Galaxy are typically tens 
of km s _1 , and may reach several hundred km s _1 , the 
ambiguity has practical relevance in the context of precise 
stellar radial velocities. 

It is seen from Eq. (J5J that the ambiguity arises when 
the quantity di b s /dt* is transformed into a velocity, i.e. 
when a model is used to interpret the data. dt b s /dt*, on 
the other hand, is a direct, model- independent relation 
between the basic events of light emission and reception. 
From an observational viewpoint, we could therefore re- 
gard the dimensionless quantity di bs/dt* as more funda- 
mental than cither v' r or v" . 

2.2. Doppler shift 

The result of a spectroscopic line-shift measurement may 
be expressed by means of the dimensionless redshift vari- 
able 

Aobs — Ao 



z 



X 



£L 



^0 — ^obs 



(6) 



2 There is an analogous problem in the definition of proper 
motion, i.e. the rate of change in direction u, but the consensus 
is that proper motion means du/dt b E , not du/dt,; cf. Sect. 14.51 



where Ao (yo) is the rest- frame wavelength (frequency). 
The redshift is often converted to a conventional velocity 
using some standard formula, the simplest being 



m A bs — Ao 

v). ' — c = cz . 

Ao 



(7) 



However, an alternative conversion is obtained by consid- 
ering the relative shift in frequency rather than in wave- 
length, viz.: 



12) V ~ v ohs 

vi > = c 



— • (8) 

v 1 + z W 

This last expres sion has tradi tionally been used in radio 
astronomy (e.g. IWalkerl fl987(l . although the practice is 
discourage d due to the risk of confusion with the earlier 
expression l|Contopoulos fc Jappe]|ll974t iMiiller fc Jappel 
Il977l) . wf 1 ^ is sometimes called the 'optical velocity' and 
v r 2) the 'radio velocity' l|Greisen et al.ll2003j) . 

For use with large ve locities, the following formula is 
often recommended (e.e. lLan3ll974|) : 



= c 



{l + zf 



1 



(1 + zf + l- (9) 

The expression is derived from the special-relativistic 
Doppler formula by assuming purely radial motion (cf. 
below) . 

Thus we have at least three different conventions for 
transforming z into v r that are more or less 'standard' in 
astronomy. From the series expansions 

(i) 

Vr = CZ 



,(3) 



c (z — z 2 + z 3 



+ •••) 



\z 2 + iz 4 



(10) 



- c [z ., . j 

it is seen that the differences between the three conven- 
tions are of second order (cz 2 ~ v 2 /c). 

Equation © ignores the star's transverse velocity (vt). 
The c omplete Do ppler formula from special relativity 
reads (|Langlll974h : 



1 + z 



1 + v r /c 



(11) 



y/1 -V 2 /C 2 

where v = (v 2 + v 2 ) 1 / 2 is the total velocity relative the 
observer. Solving for v r and expanding in powers of z and 
u = Vt/c we obtain a fourth expression for v r : 

=c[z- \{z 2 + u 2 ) + \zu 2 + l(z 4 -u 4 ) -■■■]. (12) 

Comparing with the third variant of Eq. (|10|l it is seen that 
the transverse Doppler effect is ~ v 2 /2c, i.e. typically of 
a similar size as the differences among the expressions in 
Eq. JTUJl. 

Thus, various conventions exist for converting the ob- 
served Doppler shift into a radial velocity; the differences 
are of order v 2 /c, exceeding 1 m s _1 for normal galac- 
tic stellar velocities and 100 m s _1 for high- velocity stars. 
Within the framework of special relativity (thus ignor- 
ing the many other effects discussed below), a 'rigorous' 
transformation from z to v r is possible, but only if the 
transverse velocity of the star is also known. 
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2.3. Astrometric determination or radial motion 



3.1. Gravitational redshifts 



Astrometry specialises in measuring the directions to ob- 
jects, and in particular the directional changes caused by 
the motions of the object (proper motion) and observer 
(parallax). Although such measurements primarily yield 
the distances and transverse velocities of the objects, they 
are in principle sensitive also to the radial motion of the 
objects through various second-order effects. Although the 
principles have been known for a long time, it is only with 
the high accuracies realised with space techniques that 
astrometry has become a practical possibility for radial- 
velocity measurements. 

Different methods exist for t he determination of astro- 
metric stellar radial velocities (|Dravins et alJll999^ . The 
most direct method, measuring the rate of parallax change 
as a star approaches or recedes, is still beyond realised ac- 
curacies (e.g., even for the nearby high- velocity Barnard's 
star, the expected parallax rate is only 34 //as yr _1 ), al- 
though it is expected to become measurable in the foresee- 
able future. Another method utilises that a star's proper 
motion changes as a result of its changing distance from 
the Sun ('perspective acceleration'). By combining high- 
accuracy proper motions with measurements of stellar 
positions at different epochs, radial-velocity values have 
so far been determined for some 20 stars, though only 
with modest accura cies (typically a few tens of km s _1 ; 
iDravins et alJllQQflh . However, a third method, applicable 
to nearby moving clusters such as the Hyades, whose stars 
share the same (average) velocity vector, already permits 
accuracies on the sub-km s -1 level, and for several classes 
of stars has yielded better data than have been possible 
to achieve spectroscopically. Here, parallaxes and proper 
motions are combined to determine the apparent secular 
expansion or contraction of the angle subtended by a clus- 
ter, as it is approaching or receding. Using data from the 
Hipparcos satellite mission, more than a thousand st ars 
have already been thus studied l(Madsen et alJEof)2l . a 
number to be increased when data from future astrome- 
try missions become available. 

These methods are all based on the same general prin- 
ciple: let 9 be the angular size of an object (the Earth's 
orbit as seen from the star, the distance travelled by the 
star in a given time, or the size of a stellar cluster) and 
r its distance; then the assumption is that r6 = con- 
stant, from which r = —r8/9. While the principle is simple 
enough, the question still remains whether the derivative 
( ' ) should be taken with respect to the time of observa- 
tion, or the time of light emission. Thus, also the concept 
of astrometric radial velocities needs a more precise defi- 
nition. 



3. Limitations of spectroscopic radial-velocity 
measurements 

In this Section we highlight the various issues that may 
limit the achievable accuracy in stellar radial velocities, as 
deduced from spectroscopic observations. 



The gravitational potential induced by a star's mass 
causes redshifts of all photons leaving its vicinity. Across 
the Hertzsprung-Russell diagram, the gravitational red- 
shifts change by three orders of magnitude between white 
dwarfs (some 30 km s _1 ) and supergiants (some 30 m s _1 ). 
This gravitational redshift u grav = GM/rc diminishes 
with distance from the stellar centre as r -1 . For the Sun, 
the value is 636.5 m s^ 1 for light escaping from the solar 
photosphere (r = i? ) to infinity, 3 and 633.5 m s _1 for 
light intercepted at the Earth's mean distance from the 
Sun (r — 21hRo). A solar spectral line instead formed at 
chromosphcric heights (30 Mm, say; r = 1.047?q) will have 
this redshift decreased by some 20 m s" 1 , and a coronal 
line by perhaps 100 m s _1 . 

For other stars, the shift scales as (M/M@)(R/ Rq)' 1 , 
or as (g/g & )(R/R ( ? ) ), where g is the surface gravity. Since 
R and M c an rarely be estim ated to better than ~ 5% for 
single stars l)Andersenlll99l|) . while spectroscopi c determi- 
natio ns of log g have much larger uncertainties l|Lebretonl 
2000). it is normally not possible to compute the gravi- 
tational redshift to better than 50 m s _1 for individual 
single stars. 

3.2. Effects inside stellar atmospheres 

It is well known that photospheric absorption lines in the 
solar spectrum are blueshifted by about 400 m s _1 (after 
correction for the known gravitational redshift) as a result 
of convective motions in the solar atmosphere. In the pho- 
tospheric granulation, hot (bright) and rising (blueshifted) 
convective elements contribute more photons than the 
cooler (darker) and sinking (redshifted) g as, thus caus- 
ing a net blueshift of the absorption lines (|Dravinslll982t 
lAllende Prieto k. Garcia L6pezHl998j) . Detailed modelling 
of stellar atmospheres involving 3-dimensional and time- 
dependent hydrodynamics is capable of producing syn- 
thetic spectral lines whose intensity profiles and patterns 
of wavelength displacements closel y match observations 
for at least solar-type st ars (e.g. lAsnlund et alJ [2000; 
lAllende Prieto et al.l 120021. For main-sequence stars, the 
predicted convective blueshifts range between approxi- 
mately 1000 m s^ 1 for F-type stars and 200 m s _1 for 
cooler K-types. 

The shift is however not the same for all the lines in 
a spectrum. The precise amount of shift depends on the 
strength of the absorption line (and hence on the stellar 
metallicity) , since different lines are formed at different 
atmospheric depths and thus experience different granula- 
tion contrasts. For the Sun such (observed and modelled) 
differential shifts between differently strong lines in the 
visual amount to some 200 m s _1 , but reach 1000 m s -1 
for the hotter F-type star Procyon ijAllende Prieto et alJ 



3 Using R Q = ( 6.95508 ±0.00026) x 10 8 m for the solar pho- 
tospheric radius ferown fc Christensen-Dalsgaardl [1 998) and 
GM Q = 1.327124 x 10 20 m 3 s" r 
GM Q /Rc = 636.486 ± 0.024 m s -1 . 



JStandishl Il995ft we get 
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12002ft . The shifts further depend on excitation potential 
and ionisation level (due to different conditions of line for- 
mation), and on wavelength region (due to varying gran- 
ulation contrast). Actually, in some wavelength regions, 
where the lines may originate in atmospheric layers char- 
acterised by convective overshoot (with an inverted ve- 
locity/temperature correlation), the lineshifts may change 
sign to become convective redshifts. For the Sun this can 
be observed in cores of very str ong lines in th e optical or 
generally in the far ultraviolet l|Samainlll99lll . 

Line profiles are also asymmetric, making the de- 
termination of accurate lineshifts a matter of conven- 
tion - where in the line should the shift be measured? 
Asymmetries are caused not only by convective motions 
on the stellar surface but also by asymmetric emission 
and/or absorption components (e.g., due to chromospheric 
emission or stellar- wind absorption), by microscopic pro- 
cesses causing wavelength shifts on the atomic level (e.g., 
pressure shifts), or macroscopic circumstances (e.g., grav- 
itational redshift). Further complications enter for pulsat- 
ing stars, those with expanding atmospheres, or such with 
deviant isotopic compositions. 

Since many physical effec ts thus contrib ute to the ob- 
served wavelength shifts fe.g. lDravinsll20fl.'^ . it is not pos- 
sible to deduce an accurate centre-of-mass motion simply 
from the observed differences between wavelengths in the 
source and those measured in the laboratory. 

3.3. The role of standard stars 



son is simply the physical nature of stellar spectra and 
the practical impossibility to obtain noise-free measure- 
ments. Apart from the physical effects of stellar surface 
convection and gravitational redshifts mentioned above, 
the wavelengths of stellar spectral lines depend on, La., 
factors such as the stellar rotation rate, the angle under 
which the stellar rotation axis is observed, the phase in 
a possible magnetic stellar activity cycle, and the spec- 
tral resolution and instrumental profile of the observer's 
instrument. Here, we give examples of such effects that 
are likely to limit the ultimate achievable precision for 
radial-velocity standard-stars to levels not much better 
than 0.5 km s _1 : 



3.3.1. Effects of stellar rotation 



The influence of stellar rotation has been realised, espe- 
cially for_e^jfier^ty_D_e_s_tarew 

tion llAn cjersen Nordg tronJIlQsll IVerschueren fc Davidl 
Il999t iGriffin et alJl200(t . The effects caused by the mis- 
match between a spectrum template for a slow-rotation 
star and a rapidly rotating A-typ e star may exceed 
1 km s _1 ijVerschueren fc Davidll999]l . The origin of these 
effects is the rotational line-broadening and the ensuing 
blending of spectral lines, significant when a given star is 
observed equator-on, but disappearing when viewed pole- 



Radial- velocity standards, with a supposedly known 'true' 
velocity, have long been used as objects against which to 
calibrate observations of other stars. Indeed, a goal of the 
IAU Working Group on Radial Velocity Standard Stars 
has been 'to provide a list of such standard stars whose 
velocities are known with an accuracy of 100 m s -1 .' 
Overviews of t he wo r k by this an d othe r group s are i n 
IStefanik et all Jl999l) : llldrv et alJ Jl999l) : iFekell JjjjjjjJ, 
and in various triennial reports from IAU Commission 30 
on r adial velocities , published in the IAU Transa ctions 
(e.g. lBergerorJH99i lAndersenlll999T: iRickmanlEool . 

In practical application, however, a number of depen- 
dences on the 0.5 km s _1 level have been found, in particu- 
lar on the stellar colour index, differences among different 
radial- velocity instruments, between different spectrum 
correlation masks applied on the same stellar spectrum, 
etc. The best agreement is normally found for stars of 
spectral type close to that of the Sun - naturally so, since 
the instruments and data reductions are normally cali- 
brated against the solar spectrum e.g. as reflected off mi- 
nor planets, whose motions are accurately known through 
other methods, and the procedures set up to produce con- 
sistent results at least for such a solar spectrum. 

However, not even very elaborate calibrations are likely 
to produce any 'true' standards to a much better preci- 
sion than 0.5 (or, perhaps, 0.3) km s _1 , unless a detailed 
physical model of the observed star is developed. The rea- 



Even very modest rotational velocities in sharp-lined 
late-type stars may cause significant wavelength displace- 
ments of the spectral-line bottoms and other parts of 
the line profiles used for radial-velocity determinations. 
Naively, one might expect that increased stellar rotation 
would merely smear out the line profiles and perhaps 
straighten out the bisectors which describe the line asym- 
metry. Actually, for more rapid rotation, when the asym- 
metric line components originating near the stellar limbs 
begin to affect the wings of the profile integrated over the 
stellar disk, the line asymmetries may become enhanced . 
This p henomenon w as suggested bv lGrav fe Tonerl l|l985h 
and bv lGravl (l986) and studi ed in more detail f or a simu- 
lated rapidly rotating Sun bv lSmith etaP l(l987t) . Detailed 
line profile calculations from hydro dynamic model atmo- 
spher es for stars of different types ^Dravins fc Nordlundl 
1990) show that this effect may easily displace photo- 
spheric line-bottoms by a few hundred m s _1 already in 
solar-type stars rotating with Tasini less than 10 km s . 
Also, stars may rotate not as rigid bodies, and rotation 
may be differential with respect to stellar latitude, atmo- 
spheric height, or between magnetic and non-magnetic el- 
ements. The existence of differentia l rotation is suggested 
both from studies of starspots (e.g . JWeber fc Strassmeierl 
120011 iGollier Cameron et al.ll2002r . and from analyses of 
line profiles ijReiners et alJ l2001L 
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3.3.2. Effects of stellar activity cycles 

On levels of perhaps 10-100 m s _1 , at least cooler stars 
undergo apparent radial-velocity variations during a stel- 
lar activity cycle, when differently large fractions of the 
stellar surface are covered by active regions with mag - 
netically 'disturbed' granulation (e.g.. lOrav et al]fl99ffl . 
Magnetic flux that becomes entangled among the convec- 
tive features limits the sizes and the temperature and ve- 
locity amplitudes to which these feature s develop. For so- 
lar observations, seelSnruit et al. I lll990l their Fig 1) an d 
ISchmidt et all <|l98S|) : for theory. see lBercik et all <|l99S|) . 
The resulting spectral line asymmetries are changed in 
the sense of smaller asymmetries and smaller wavelength 
shifts in the active regions (e .g. Ilmmerschitt fc Schroterl 
Il989[ [Brandt fe Solanki||l i)90l and references therein). 

~ ^vhig^ton"e^dTi^99Ti followed the full-disk asymme- 
tries of Fei lines during a full 11-year solar activity cy- 
cle, finding cyclic variations in the line asymmetry with 
an amplitude of about 20 m s ; presumably the corre- 
sponding absolute shifts are (at least) of a similar size. 
Indeed, variations on this order can be predicted from 
spatially resolved observations of line profiles in active re- 
gions, weighted with the cyclically changing area coverage 
of active regions during an activity cycle. The effects can 
be expected to increase (to perhaps 50 m s _1 ) for younger 
and chromosphcrically m ore active stars, e.g., F - and G- 
type ones in the Hyades l)Saar fc Donahuelll997l) . 

Since the amount of convective lineshift differs among 
different types of spectral lines and between different spec- 
tral regions, also the activity- induced changes in this shift 
must be expected to differ. While the identification of such 
differences could be important to find lines whose sensi- 
tivity to stellar activity is smaller (thus being better diag- 
nostics for exoplanet signatures) or greater (being better 
diagnostics for magnetic activity), such data are not yet 
available (and may indeed require a stringent definition 
of the radial- velocity measure to permit intercomparisons 
between observations at different epochs). 

Besides these cyclic changes, there are shorter-term 
fluctuations (on a level of perhaps 20-30 m s _1 ) in the ap- 
parent radial velocity of stars, which often are greater in 
stars with enhanced chromospheric activity. Presumably, 
this reflects the ev olution and changing area coverages of 
^active^regions (e.gj. ISaar et aDll998USaar fe Fischerll200fj 

3.3.3. Effects caused by starspots and surface 
inhomogeneities 

Greater effects are present in spotted stars whose pho- 
tometric variability indicates the presence of dark spots 
across the stellar surface. The amplitude of variations 
expected from photometrically dark spots is on the or- 
der of 5 m s -1 for solar-age G-type stars, increasing to 
perhaps 30-50 m s" 1 for younger and more active stars 
1 Saar fc Donahue! Il997t iHatzes! Eool . In more heavily 
spotted stars (such already classically classified as vari- 



ables) , the technique of Doppler imaging exploits the vari- 
ability of spectr al-line profiles to rec onstruct stellar sur- 
face maps (e.g.. IPiskunov et"aflll990h . but obviously any 
more accurate deduction of the radial motion of the stel- 
lar centre-of-mass from the distorted spectral lines is not 
a straightforward task. 

Surface inhomogeneities causing such line distortions 
need not be connected to photometrically dark (or bright) 
spots, but could be chemical inhomogeneities across 
the stellar surface (with locally different line-strengths) 
or just patches of g ranulation with different structure 
!TonerfcGra,vlll98ii. 

3.3.4. Effects caused by the finite number of granules 

Even the spectrum of a hypothetical spot-less and non- 
rotating star with precisely known physical and chemi- 
cal properties will probably still not be sufficiently stable 
to serve as a 'standard' on our intended levels of accu- 
racy. One reason is the finite number of convective fea- 
tures (granules) across the stellar surface. For the Sun, a 
granule diameter is on the order of 1000 km, and there 
exist, at any one time, on the order of 10 6 such granules 
on the visible solar disk. The spectrum of integrated sun- 
light is made up as the sum of all these contributions: to 
make an order-of-magnitudc estimate, we note that each 
granule has a typical velocity amplitude of 1-2 km s -1 . 
Assuming that they all evolve at random, the apparent 
velocity amplitude in the average will be this number di- 
vided by the square root of 10 6 , or 1-2 m s _1 . This 'astro- 
physical noise' caused by the finite number of granules is 
a quantity that is becoming measurable also in solar-type 
stars in the form of an excess of the power spectrum of 
spectral-line variability at temporal frequencies of some 
mHz, corresp onding to granular li fetimes on the order of 
ten minutes jKieldsen et alJll999l). Although it can als o 
be modelled theoretically (e.g. iTrampedach et alJll99Sfl . 
such modelling can only predict the power spectrum and 
other statistical properties, not the instantaneous state of 
any star. 

The number of granules across the surfaces of stars 
of other spectral types may be significantly smaller than 
for the Sun, and the resulting 'random' variability corre- 
spondingly higher. For supergiants, it has been suggested 
than only a very small number of convective elements (per- 
haps only a few tens) coexist at any given time, but even 
a star with 10 4 granules would show radial-velocity flick- 
ering an order of magnitude greater than the Sun. 

The effects are qualitatively similar for other types 
of statistically stable structures across stellar surfaces, 
e.g. p-mode oscillations, where various surface regions 
on the star are moving with varying vertical velocities. 
Since their averaging across the star does not fully can- 
cel, they are detectabl e as lineshift variations on a level 
of several m s" 1 (e.g . lHatzesI IT9f)3 iBedding et aT1l200lt 
iFrandsen et alJEf)02lf . demonstrating another application 
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of precise radial- velocity measurements, as well as the lim- 
itations to stellar wavelength stability. 



opacities make the deeper layers visible, with perhaps 
greater convective amplitudes. 



3.3.5. Effects caused by exoplanets 

Intrinsically stable stars may show variability on the 10- 
100 m s _1 level induced by orbiting exoplanets. Of course, 
this is a true radial-velocity variation, but it practically 
limits the selection of such stars as radial-velocity 'stan- 
dards', since their use on the m s _1 level would require 
detailed ephemerides for their various exoplanets. It can 
be noted that 51 Peg used to be a radial- velocity standard 
star! 



3.3.6. Effects of instrumental resolution 

A different type of wavelength displacements is introduced 
by the observing apparatus, in effect convolving the pris- 
tine stellar spectrum with the spectrometer instrumental 
profile. Since all stellar spectral lines are asymmetric to 
some extent, their convolution with even a perfectly sym- 
metric instrumental profile of an ideal instrument pro- 
duces a different asymmetry, and a different wavelength 
position e.g. of the line-bottoms. Quantitative calculations 
demonstrate how such effects reach 50 m s -1 and more, 
even for high-resolution instrum ents l)Brav fc Louehheadl 
Il978l iDravins fc Nordmndlll990h . Although this is 'only' 
a practical limitation which, in principle, could be cor- 
rected for if full information of the instrumental response 
were available, this is not likely to be possible in practice. 
For example, differences of many tens of m s _1 in mea- 
sured lineshifts may result between spectrometers with 
identical spectral resolutions (measured as full width at 
half-maximum), but whic h differ only in their amounts 
of diffuse scattered light l)Dravinslll987|) . Of course, even 
greater lineshifts could be caused by asymmetric instru- 
mental prof iles. Inst r umen tal effects in spectroscopy are 
reviewed bv lDravinsI l|l994l) while methods for calibrating 
instru ment al profiles are disc ussed by, e.g., IValenti et alJ 
l|l995(l and lEndl etall l|2000j) . 



3.3.7. Effects of instrumental design 

For observational modes not involving analyses of highly 
resolved line profiles, but rather statistical functions such 
as the cross-correlation between spectral templates, a se- 
ries of other instrumental effects may intermix with stel- 
lar ones. For example, if a detector/template combination 
predominantly measures a signal from the blue spectral 
region, it may be expected to record a somewhat greater 
spectral blueshift in cool stars since convective blueshifts 
generally increase at shorter wavelengths (where a given 
temperature contrast in surface convection causes a rela- 
tively greater brightness contrast). A red-sensitive system 
may give the opposite bias, unless it is sensitive into the 
infrared, where the generally smaller stellar atmospheric 



3.3.8. Conclusions about spectroscopic 'standard' stars 

Physical and instrumental effects, such as those listed 
above (and others, such as errors in laboratory wave- 
lengths), imply that there most probably do not exist any 
stars whose spectral features could serve as a real stan- 
dard on precision levels better than perhaps 300 m s . Of 
course, for poorer precisions - perhaps around 0.5 km s _1 
- various standard sources, including the solar spectrum, 
may continue to be used as before. However, in order to 
deduce true velocities to high accuracy, all spectral obser- 
vations - of 'standard' stars and others - must undergo 
a detailed physical modelling of their emitted spectrum, 
and of its recording process. 



3.3.9. Possible future astrometric standard stars 

The recently realised accurate determination of stellar ra- 
dial motions through astrometric measurements opens the 
possibility of having also radial-velocity standards deter- 
mined independent of spectroscopy. While the ultimate 
limitations in thus obtainable accuracies have not yet been 
explored (e.g., what is measured in astrometry is the pho- 
tocentre of the normally unresolved stellar disk, whose co- 
ordinates may be displaced by starspots or other features) 
there appear to be no known effects that in principle would 
hinder such measurements to better than 100 m s _1 (or 
even 10 ms -1 ), at least for some nearby stars. This will re- 
quire astrometric accuracies on the microarcsecond level, 
and possibly extended periods of observations, but these 
are expec ted to be reachable in future space astrometry 
missions (|Dravins et alJll999j) . 

3.4. Conclusions 

From the above discussion it is clear that numerous effects 
may influence the precise amount of spectral- line displace- 
ments. Among these, only local effects near the observer 
(i.e., within the [inner] solar system) can be reliably cal- 
culated and compensated for. In particular, these depend 
on the motion and gravitational potential of the observer 
relative to the desired reference frame, normally the solar- 
system bary centre. 

Therefore, spectroscopic methods will not be able, in 
any foreseeable future, to provide values of stellar radial 
motion with 'absolute' accuracies even approaching our 
aim of 1 m s . Of course, this does not preclude that 
measurement precisions (in the sense of their reproducibil- 
ity) may be much better and permit the detection of very 
small variations in the radial velocity of an object (whose 
exact amount of physical motion will remain unknown) 
in the course of searching for, e.g., stellar oscillations or 
orbiting exoplanets. In order to enable further progress in 
the many fields of radial- velocity studies, it seems however 
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that more stringent accuracy targets have to be defined, 
so that future observational and theoretical studies have 
clear goals to aim at. 

Astrometric radial velocities do not appear to have the 
same types of limitations as those deduced from spectro- 
scopic shifts, and more lend themselves to definitions that 
can be transformed to 'absolute' physical velocities. These, 
however, must be stringently defined since different plau- 
sible definitions differ by much more than our desired ac- 
curacies. 



4. Kinematic and astrometric radial velocity 

We will now more thoroughly scrutinise the various ge- 
ometric effects entering the concept of 'radial velocity', 
aiming at definitions that are consistent at an accuracy 
level of 1 m s _1 . This requires first that a system of tem- 
poral and spatial coordinates is adopted; then that the 
relevant parts of the astronomical event (Fig. 1) are mod- 
elled in this system, consistent with general relativity at 
the appropriate accuracy level; and finally that suitable 
conventions are proposed for the parameterisation of the 
event. 

The modelling of astrometric observations within a 
general-relati vistic framewo r k has been t reat ed in text- 
books such as lMurravl l|l983|) . ISoffe]| l|l989|) and lBrumbe73 



several papers (e.g. Stumnfl 


1985; Backer & ] 


ffellin ff sJl98f4 


iKlioner & Koneikinl Il992l 


iKlionerl 12000b 




2(U(). Much 



of the mathematical development in this and the next 
Section is directly based on these treatments, but adapted 
in order to present a coherent background for the defini- 
tion and explanation of the radial- velocity concepts. 

4.1. Coordinate system (BCRS) 

Subsequently, t and r = (x, y, z) denote coordinates in the 
Bary centric Celestial Reference Syst em (BCRS) ado pted 
by the IAU 24th General Assembly jRickmanll200ll) and 
discussed by iBrumberg fc Crotenl l|200lf) . The temporal 
coordinate t is known as the Barycentric Coordinate 
Time (TCB). The BCRS is a well-defined relativistic 4- 
dimcnsional coordinate system suitable for accurate mod- 
elling of motions and events within the solar system. 
However, it also serves as a quasi-Euclidean reference 
frame for the motions of nearby stars and of more distant 
objects, thanks to some useful properties: it is asymptot- 
ically flat (Euclidean) at great distances from the Sun; 
the directions of its axes are fixed with respect to very 
distant extragalactic objects; and the origin at the solar- 
system barycentre provides a local frame in which nearby 
(single) stars appear to be non-accelerated, as they experi- 
ence practically the same galactic acceleration as the solar 
system. The axes are aligned with the celestial system of 
right ascension and declination as realis ed, for exam ple, 
by the Hipparcos and Tycho Catalogues ljESAlll997l) . 



4.2. The light-time equation 

As emphasised in Sect. l2.ll the relation between the events 
of light emission and light reception is fundamental for de- 
scribing the astronomical event resulting in a geometric or 
spectroscopic observation. The two events are connected 
by the light-time equation, from which the required spatial 
and temporal transformations may be derived. 

In the BCRS, let r»(i) describe the motion of the star 
and r Q bs(i) that of the observer. Assume that a light signal 
is emitted from the star at time , when the star is at the 
spatial coordinate r» = ?%(£») and its coordinate velocity 
is v± = (dr*/dt) t =t t - Assume, furthermore, that the same 
light signal is received by the observer at time i bs, when 
the observer is at the spatial coordinate r b s = »"obs(iobs) 
and its coordinate velocity is v b s = (dr b s /d£)t = t oba . The 
light-time equation can now be written 



^obs ^* C ^obs | ~\~ T (^* 5 > ^obs ) ^*obs) 



(13) 



Here, c = 299 792 458 m s _1 is the conventional speed of 
light, | | denotes the usual (Euclidean) vector norm, and 
7 is the relativistic delay of the signal along the path 
of propagation from star to observer. The delay term is 
required to take into account that the coordinate speed 
of light in the presence of a gravitational field is less 
than c in the adopted metric, so that the first term in 
Eq. I|13fl gives too small a value for the light travel time. 
For present purposes it is sufficient to describe the grav- 
itational field by means of the total Newtonian poten- 
tial $(r) relative to the BCRS. For instance, the gravita- 
tional field of the solar system is adequately described by 
$(r) = GJ^i Mi\f — T*i | — 1 , where G is the gravitational 
constant, and the sum is taken over the different solar- 
system bodies having (point) masses Mi located at coordi- 
nates To first order in c~ 2 , the coordinate speed of light 
in the BCRS metric is given by |dr/dt| = c(l-2$/c 2 ). The 
time delay per unit length is therefore 2<3>/c 3 . Integrating 
this quantity along the light path (which for this calcu- 
lation can be taken to be a straight line in the BCRS 
coordinates) gives 4 
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E 



2GM, 



hi 



k (v* Vi 



r* - n 



(14) 



fc'(r-obs - ri) + |r obs - r 
k is the coordinate direction from the observer to the star, 
k = |r» - robsp 1 ( r * - fobs) • (15) 

Subsequently, we shall mainly consider the gravita- 
tional field of the Sun (index = 0), for which 2GM Q c~ 3 ~ 
9.85 /is. For objects as distant as the stars we can neglect 
Tq compared with r* in the numerator of the argument to 
the logarithm in Eq. Ijl4jl . Moreover, k is practically par- 
allel with r„, so the numerator becomes simply 2r„, where 
r* = \r*\. (The error introduced by this approximation is 

4 Since the perturbing bodies move during the light propa- 
gation, r, should be taken to be the position at the time of 
closest approach of the photon to the perturbing body. For a 
rigorous treatment, see lKopeikin fc SchafeJ lll999ll . 
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< 10~ 16 s for r* > 1 pc.) The denominator varies depend- 
ing on the relative positions of the Sun, observer and star, 
but is typically of the order of the astronomical unit (A) 
for an observer on the Earth. Thus 7 ~ (10 fis) ln(2r*/A), 
or ~ 100 /zs for the nearest stars, ~ 200 [is at i\ = 10 kpc, 
and ~ 300 /is for objects at cosmological distances. This 
slowly varying delay of a few hundred microseconds suf- 
fered by the light while propagating from the star to the 
solar system is generally ignored. Indeed, it would hardly 
make sense to try to evaluate it, since the gravitational 
delays caused by other bodies (in particular by the star 
itself) are not included. 

However, there is also a rapidly varying part of 7 in 
Eq. 1(14(1 , caused by the motion of the observer with respect 
to the Sun. In order to separate the rapidly varying part 
of the delay from the (uninteresting) long-range delay, we 
write 



E2GMi ( k'(r ohs - n) + \r ohs - r % 



.4 



where 



7* 



2G 



M, 



In 



2r» 



(16) 



(17) 



is practically a constant for the star (cf. below). Clearly 
any constant length could have served instead of A to 
separate the terms in Eq. 1)16(1 . Using the astronomical 
unit for this purpose is just an arbitrary convention. 

From the light-time equation we can now determine 
the relation between the coordinate time interval dt* (e.g. 
representing one period of emitted radiation) and di b s 
(the corresponding period of coordinate time at the ob- 
server). Writing |r* — r b s | = fc'(r* — r Q bs) we have 



diobs - di* = fc'(dr% - dr obs )c 1 + d7 , 



from which 



1 + k'v ohs c 1 - dj/dt 



obs 



obs 



1 + k'v* 



(18) 



(19) 



From Eqs. I(16|) - ((17(1 it follows that d7/di b s is the sum 
of two terms, the first of which ~ (10 /us)p/r*, where 
p = drv/d^B is the astrometric radial velocity of the star 
defined below (Sect. 14.5(1 . For |p| < c and r, > 1 pc this 
term is < 10~ 13 and therefore negligible. The second term 
can be evaluated e.g. for an observer in circular orbit (at 
1 AU) around the Sun. It reaches a maximum value of 
< 10~ 9 when the observer is behind the Sun so that the 
light ray from the object just grazes the solar limb. The 
simple expression 



di* 



obs 



1 + k'v ohs c 1 

1 + k'v* c- 1 



(20) 



is therefore always good enough to a relative accuracy 
better than 10~ 9 (< 0.3 m s _1 in velocity). 



4.3. Barycentric time of light arrival 

Since distances to objects beyond the solar system are gen- 
erally not known very accurately, it would be inconvenient 
to use the time coordinate for describing observations 
of phenomena that occur at such great distances. Instead, 
it is customary to relate the observed events to the time 
scale of the observer. For very accurate timing, e.g. as re- 
quired in pulsar observations, one must take into account 
both the geometrical (R0mer) delay associated with the 
observer's motion around the solar-system barycentre, and 
the rclativistic (Shapiro) delay caused by the gravitational 
field of bodies in the solar system. 

We define the barycentric time of light arrival as 



tB = t* + r* c 1 + 7* 



(21) 



where 7* is given by Eq. 1(17(1 . That is, ie is the time of 
light emission delayed by the nominal propagation time 
to the barycentre (r*/c) plus that part of the relativistic 
delay which is independent of the observer. By means of 
Eqs. ((14(1 and 1(16(1 we find that the barycentric time of 
light arrival becomes 



tabs + (r* - R)c 



2GM; 



111 



k'{r ohs 



^obs 



.4 



(22) 



where R = r* — r b s | is the topocentric coordinate dis- 
tance to the star. With r* ~ ttr», where u is the barycen- 
tric coordinate direction to the star, the R0mer delay can 
be expanded to give 



(r* - R) c 



u'r ohs \u x r obs | 2 , (u'r ohs ) 



2cr* 



2cr? 



(23) 



For r, > 1 pc and r b s ~ 1 AU the maximum amplitudes 
of the three terms are ~ 500 s, 1 ms, and 4 ns, respectively; 
neglected terms are of order < 10~ 13 s. 

While Eq. I(21|) formally defines the barycentric time 
of light arrival, it is clear that Eqs. I(22l) - ((23(l must in prac- 
tice be used to calculate it for any given observation. In 
principle this also requires that the distance r» is known, 
but only to a moderate precision. In many practical situ- 
ations, the curvature terms in Eq. ((23(1 [depending on r^ 1 
and r~ 2 ] can be neglected. 

4.4. Definition of kinematic parameters 

Within the BCRS we define the kinematic parameters 
of a star to be its coordinate at time £», and the 
instantaneous coordinate velocity at the same instant, 
u* = dr*/dt\t=t, ■ In a time interval around t* we have 



r*(t) = r*(t*) + (t - t*)v* + 0(t - i*) 2 



(24) 



The six components of the phase-space vector (r„, v*) are 
the relevant elements for studies of galactic kinematics and 
dynamics, for example integration of galactic orbits (after 
transformation to a galactocentric system). 
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We define the kinematic radial velocity as the compo- 



nent of i>* along the barycentric direction u = r* 



(25) 



The perpendicular component of the coordinate velocity 
is the kinematic tangential velocity, v t — — uv r . The 
kinematic radial and tangential velocities are equivalent 
to the 'true' rad ial and tangential velocities introduced by 
iKlionerl ll2000bh . 

4.5. Definition of astrometric parameters 

The six components of (j%, v*) are not directly observable 
but can in principle be derived from astrometric obser- 
vations of the star. There is consequently an equivalent 
set of six astrometric parameters, which we now define. 
Following the discussion in Sect. 14.31 the astrometric pa- 
rameters are considered as functions of t^. Writing the 
barycentric coordinate of the star as 

r* = itr* (26) 

we define the celestial coordinates (a, 5) of the star at the 
epoch ie by means of the components of the unit vec- 
tor u in the BCRS. This gives the first two astrometric 
parameters. The third one is parallax, which we define as 



.4 



(27) 



rcf. lKlionerll2000b() . The rate of change of the barycentric 
direction is the proper-motion vector, 



A 4 



dit 



(28) 



from which the proper motion components p a * = a cos 8 
and fig = 5 follow (the dot signifies differentiation with 
respect to £b)- These five parameters are practically iden- 
tical to the standard astrometric parameters used, for in- 
stance, in the Hipparcos Catalogue. 5 The sixth astromet- 
ric parameter is the rate of change in barycentric coordi- 
nate distance, which we call the astrometric radial veloc- 
ity: 



dr* 
die" 



(29) 



The term is motivated because of the exact analogy 
with the definition of the (astrometric) proper motion in 
Eq. I|28p. Observationally, the astrometric radial velocity 
can in principle be determined e.g. from the secular change 

There is however a subtle difference, in that proper motions 
in the Hipparcos Catalogue are formally defined as du/dTs, 
where Tb is the barycentric time of light arrival expressed on 
the Terrestrial Time (TT) scale. The TT scale is essentially the 
observer's proper time, and differs from the coordinate time 
(TCB) used in Eg. bv the averag e factor (AT /At) ~ 1 - 
1.55 x 10" 8 Jlrwin fc Fukushimall 999l [cf. Eqs. |£U and gSJ]. 
The Hipparcos proper motions should therefore be multiplied 
by 0.9999999845 in order to agree with the present definition. 



in parallax, p = A(A/ir)/At B l|Dravins et alJll999|) . The 
vector /zr* may be called the astrometric tangential ve- 
locity. The astrometric radial and tangential velocities are 
equivalent to the ' apparen t ' radia l and tangential veloci- 
ties introduced bv lKlionerl ll2000bh . 

It should be noted that the astrometric radial veloc- 
ity is conceptually quite different from the spectroscopic 
radial- velocity measure to be defined in Sect. |S] The as- 
trometric radial velocity refers to the variation of the co- 
ordinates of the source, and therefore depends on the cho- 
sen coordinate system and time scale. By contrast, the 
outcome of a spectroscopic observation is a directly mea- 
surable quantity and therefore independent of coordinate 
systems. 

4.6. Transformation between kinematic and 
astrometric parameters 

The barycentric coordinate of the star is immediately 
derived from the astrometric parameters by means of 
Eqs. (gSHlTj!, viz.: 

r* = u{A/n) . (30) 
Taking the derivative of Eq. i|21[) with respect to t-Q gives 



d<B c 



(31) 



The second term between parentheses is < 10 13 for r* > 
1 pc. To sufficient accuracy we have therefore 



di* 



1 



P 



die c 

If Eq. (|26|) is differentiated with respect to we find 



(32) 



dit 



dr» 



k* = — r* + it— 

dr» dr* 



dit dr* \ die 
r* + it I 

d^B dte / ch» 



= (/ir, + up) (1 - p/c)- 1 . (33) 
Separating the radial and tangential components we have 
v r = p(l- p/c)- 1 (34) 
and 

v t = p.r^ (1 - p/cy 1 . (35) 

Equations l|30|l and (|33|I " H35(I provide the complete trans- 
formation from astrometric to kinematic parameters. For 
the inverse transformation, we immediately obtain u and 
tt from r» by means of Eqs. (E1J) and (EZJ). Multiplying 
Eq. i.VMi scalarly with it gives 



P 



tt'i>* 



1 + u'v* c 1 
from which finally 

TT V* — Ull'v* 

~ A 1 + ii'd* c" 1 



v r (l + v r /c)~ 



(tt/A) v t (1 + v r /c) 



(36) 



(37) 
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The kinematic quantities u*, v r and Vt are coordinate 
speeds of the object and therefore physically bounded by 
the local coordinate speed of light, which in the BCRS 
far away from the Sun is very close to c. The astrometric 
radial velocity p and the astrometric tangential velocity 
£tr*, on the other hand, are apparent quantities which 
may numerically exceed the speed of light. This is so be- 
cause the denominator 1 + v r /c in Eqs. I|36|l and l|37(l can 
become arbitrarily small for an object moving at great 
speed towards the observer. Thus, p < — c if v r < ~\c, 
while | \xr* > c if v t > c + v r . The effect is equivalent to 
the standard kinematic explanation of the superluminal 
expansion observ ed in many extragalactic radio sources 
I Bedford et aTlll977llVerrneulen fc Cohenlllflfll. 



5. The spectroscopic parameter: Barycentric 
radial-velocity measure 

We have found that the naive notion of radial velocity 
as the line-of-sight component of the stellar velocity is 
ambiguous already in a classical (non-relativistic) formu- 
lation. In a relativistic framework the observed shift de- 
pends on additional factors, such as the transverse velocity 
and gravitational potential of the source and, ultimately, 
the cosmological redshift. Since these factors are generally 
not (accurately) known to the spectroscopic observer, it 
is impossible to convert the observed shift z into a precise 
kinematic quantity. 

What can be derived from spectroscopic radial- 
velocity measurements is the wavelength shift zb cor- 
rected for the local effects caused by the motion of the 
observer and the potential field in which the observation 
was made. 6 For convenience, this shift may be expressed in 
velocity units as czb, where c is the conventional value for 
the speed of light. Although this quantity approximately 
corresponds to radial velocity, its precise interpretation is 
model dependent and one should therefore avoid calling 
it 'radial vel ocity'. The term radial -velocity measure was 
proposed bv lLindegren et al . (1999). and accepted in the 
later IAU resolution, emphasising both its connection with 
the traditional spectroscopic method and the fact that it 
is not quite the radial velocity in the classical sense. 

5.1. The observed spectral shift 

In the previous sections the time coordinates of the various 
events were all expressed on a single time scale t, i.e. the 
Barycentric Coordinate Time (TCB). As we now move on 
to consider spectroscopic measurements, it is necessary to 
include proper time (r) in our discussion. The reason for 
this is that the atomic transitions generating spectral lines 
can be regarded as oscillators or clocks that keep local 
proper time. The measurement of spectroscopic line shifts 



6 The index b for barycentric signifies that - in contrast to 
the case in cosmology - the shift (and velocity czb) is referred 
to the solar-system barycentre, not the rest-frame defined by 
the cosmological microwave background. 



is therefore equivalent to comparing, by means of light 
signals, the apparent rates of two identical atomic clocks, 
one located at the source and the other at the observer. 

Let t„ be the proper time at the source of radiation, 
and r bs the proper time of the observer. Suppose that 
n = i^dr* cycles of radiation are emitted at frequency v* 
in the interval dr* of proper time at the source. Let us also 
suppose that the n cycles are received in the interval dr b s 
of proper time of the observer, who consequently derives 
the frequency v Q h s = n-/dr b s = zA.dT*/dTbbs- In terms of 
wavelength (A = cjv) the observed spectroscopic shift z Q ^ s 
is 

A bs _ v * _ dT bs 

A* t'obs dr* 



1 + Zobs = 



(38) 



where A* (— Xub) is the rest-frame wavelength of the spec- 
tral line. A spectroscopic lineshift measurement is there- 
fore equivalent to a direct comparison of the proper time 
scales at the source and observer. We need to relate these 
proper time scales to the coordinate time t used in previ- 
ous sections. 

The relation between proper time and coordinate time 
is defined by the adopted metric. For the Barycentric 
Celestial Reference System the accurat e transformation 
can be found for instance in lPetitl (|2000(). For the present 
applications we can ignore terms of order c -4 , leading to 
the simple transformation 



dr 
"dT 



(t,r,v) 



1 



1 



$(r) 



2 



(39) 



in which $(r) is the Newtonian potential introduced in 
Sect. |4~2*1 and v = dr/dt the coordinate velocity. This 
transformation applies both to the source (dr*/di*) and 
the observer (dr b s /di bs)- Thus, using Eq. l|2TJ|l . we have 

dr Q bs dr bs di Q bs di* 



dr* 



diobs di* dr* 



= 1- 



\Vobs\ 



x 1 



fobs 

c 2 2c 
k'v 



X 1 + 



fc'^obs 



X 1 n ~ 



2c z 



(40) 



According to Eq. l|3*5|) this equals the observed wavelength 
ratio 1 + z bs- 

5.2. Reduction to the barycentre 

The first two factors on the right-hand side of Eq. I|40|l 
depend on local conditions such as the motion of the ob- 
server in the BCRS and the gravitational potential of the 
observer. These vary between different times and locations 
of observers in the solar system, but they are also com- 
putable to high accuracy from known data, including the 
barycentric position and velocity of the observer. The last 
two factors, on the other hand, contain several quanti- 
ties that cannot be uniquely separated based on spectro- 
scopic observations. They depend on the line-of-sight com- 
ponent of the star's coordinate velocity (k'v*), but also on 
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the gravitational potential in the light-emitting region and 
(through |u*| 2 ) on the tangential velocity of the star. 

Let zb be the spectral shift corrected for the local, 
accurately computable effects, i.e. reduced to the solar- 
system barycentre. In the approximation of Eq. I|4(J|I we 
have: 



1 + Z B = (1 



^obs 



1 - 



obs 



Fobs 

2c 2 



k fobs 



(41) 



It is important to note that the unit vector k in Eq. (|41|l 
is the coordinate direction to the star given by Eq. I|15|) . 
not the observed (aberrated and refracted) direction. 

We now define barycentric radial- velocity measure as 
the quantity czb, where c = 299 792 458 m s _1 . For con- 
venience, czb is expressed in velocity units through mul- 
tiplication with the constant c. The radial-velocity mea- 
sure therefore obtains physical dimensions of SI metres 
per SI second. 7 The epoch of any spectroscopic observa- 
tion should be given as the corresponding barycentric time 
of light arrival CSect. I4.3|l . 

For an observer on the surface of the Earth (index 
= ©) we have ($ obs ) ^ GM Q /A + GM m /R m ~ 8.934 x 
10 s m 2 s -2 and (w bs) — 29785 m s _1 , so that the second 
factor on the right-hand side of Eq. (|41[) is on the average 



1 - 



$obs 



Fobs 

2c 2 



1 + 1.550 x 10" 



(42) 



In velocity units, the respective contributions from the 
solar gravitational potential, the Earth's own gravita- 
tional potential, and the Earth's velocity correspond to 
3.0, 0.2 and 1.5 m s — 1 . The main variation in this fac- 
tor comes from the annual variation in the observer's 
speed and distance from the Sun due to the eccentricity 
(e ~ 0.01671) of the Earth's orbit. The resulting ampli- 
tude is 2eGM Q /c 2 ~ 3.3 x 10~ 10 , or 0.1 m s" 1 in velocity 
units. For an Earth-bound observer, therefore, we may to 
sufficient accuracy (~ 0.1 ms _1 ) use the average factor in 
Eq. I|42l) when reducing the observed shift to the barycen- 
tre. 

5.3. Interpretation of the radial-velocity measure 

To first order, czb corresponds to the 'classical' spectro- 
scopic radial velocity. However, we emphasise that the 
radial- velocity measure is just a quantification of the spec- 
troscopic shift, not of physical velocity. Indeed, the inter- 
pretation of czb in terms of a kinematic or astrometric 
radial velocity is non-trivial and perhaps even impossible 
at the desired accuracy level. This is compounded by the 
additional effects discussed in Sect. |31 e.g. from motions 
in the stellar atmosphere, pressure shifts, and cosmolog- 
ical redshift. These effects were ignored in Eq. H40[> and 

7 Naturally, czb can be expressed in m s _1 or km s _1 ac- 
cording to convenience, and for cosmological velocities the di- 
mensionless measure zb may be preferred. 



we now introduce an extra factor 1 + X to take them into 
account. The barycentric radial- velocity measure is then 
given by 



1 + z B = 1 



k'v, 



1 - 



2c 2 



(1 + X). (43) 



Only if X, and the tangential velocity i>t are known 
to sufficient accuracy is it possible to derive k'v* from 
the barycentric radial- velocity measure. Using also the dis- 
tance information, the kinematic radial velocity v r = u'v* 
follows, and hence the astrometric radial velocity from 
Eq. (|36|l . Accurate transformation of czb to v r or p is 
therefore possible only in special circumstances. 



6. The IAU resolutions, and their application 

Based on the above discussion, and an interchange of opin- 
ions in the community during a few years, two resolutions 
for the stringent definition of spectroscopic and astro- 
metric radial-velocity concepts were adopted at the IAU 
XXI Vth General As sembly held in Manchester, August 
2000 l|Rickmadl200i| . Their full text is in the Appendix; 
in this Section we comment on their practical implications. 



6.1. 'Barycentric radial-velocity measure': A stringent 
definition for spectroscopic measurements 

Briefly, the first resolution defines the barycentric radial- 
velocity measure czb as the result of a spectroscopic mea- 
surement of line shifts; here c is the speed of light and Zb 
the wavelength shift referred to the solar-system barycen- 
tre. The definition avoids any discussion on what the 'true' 
radial velocity of the object would be. The transforma- 
tion between the spectroscopically determined barycen- 
tric radial-velocity measure czb and the physical velocity 
of the object is model-dependent and cannot be treated in 
isolation from, e.g., the tangential motion (cf. Sect. 15. 3|) . 

The definition implies that high-accuracy radial- 
velocity observations should be reduced to the solar- 
system barycentre according to procedures based on gen- 
eral relativity and using constants and ephemerides con- 
sistent with the required accuracy. 

6.2. Practical application of the spectroscopic 
definition 

For work at modest accuracies, the new definition implies 
no change of existing procedures, nor of any published 
radial-velocity values. The use of the 'barycentric radial- 
velocity measure' will only be required when absolute ac- 
curacies on the sub- km s _1 arc needed. However, its use 
permits to exploit radial-velocity measurements for new 
classes of tasks, such as studying the physical processes in 
stellar atmospheres exemplified in Sect.EI 
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6.2.1. Publishing observed wavelength shifts 

Traditionally, most published values for (stellar) radial ve- 
locities have been transformed by the observer to some 
'standard' system: instrumental; calibrated against stan- 
dard stars; against the spectrum of sunlight; or other. 
The traditionally reached precision has often been on the 
order of 1 km s _1 or perhaps slightly better. However, 
given that the recently much improved measuring preci- 
sions have begun to reach levels of m s _1 , while absolute 
calibrations are only some order of magnitude worse, this 
procedure should change. The main point in the definition 
of the 'radial-velocity measure' is that highly precise ob- 
servations should be published (also) without the observer 
trying to calibrate them against purported 'standard' ob- 
jects in an effort to deduce the objects' physical veloci- 
ties. Rather, the observations should be reduced to the 
solar-system barycentre, as detailed in the IAU resolution, 
and any subsequent interpretation of these observed wave- 
length or frequency displacements in terms of the object's 
motion, or other effects, should be made separately. The 
uncertainties in any attempted deduction of the physical 
velocity are likely to be much greater than those currently 
reachable in measurements of the wavelength shifts them- 
selves. Therefore, any precise observational data are likely 
to become corrupted by applying such model-dependent 
'corrections', rendering the data useless for possible more 
sophisticated analyses in the future. 

The 'barycentric radial-velocity measure' is a quan- 
tity that may be quite different for different spectral 
lines in the same star, or for different portions of the 
sa me spectral feature. Nice e xamples of this are seen 
in lAllende Prieto et alJ (j2Q02), where weak absorption 
lines in the spectrum of Procyon are observed to be sys- 
tematically blueshifted by almost 1000 m s -1 from the 
strong lines. Further, the wavelength positions of the line- 
bottoms are blueshifted by some 200 m s _1 relative to 
those of the line flanks closer to the continuum. These 
particular effects can be well modelled by hydrodynamical 
model atmospheres, and are found to be caused by correla- 
tions between temperature and vertical velocity in stellar 
surface convection. Thus, precise radial-velocity measures 
may be used as a nove l tool to diagnose stellar hydrody- 
namics l|Dravinsl f2003'l. provided the data have not been 
corrupted by futile attempts to 'calibrate' the apparent 
velocities. 

High-precision spectrometers often use some spectral 
template with which the observed spectrum is cross cor- 
related in order to obtain a wavelength shift. For any one 
stellar spectrum, the resulting wavelength shift will nat- 
urally depend on the exact properties of each different 
template (which portions of what types of lines are being 
selected), in which particular wavelength region are the 
measurements being made, as well as on other parameters 
llGriffin et alJl200fj IVerschueren fe Davidlll999t iGullberd 
19991 lOullberg fc Lindegrenll2002h . To retain the maxi- 
mum amount of information and permit later physical 
modelling, the barycentric radial-velocity measure should 



be given together with details of the spectral template and 
the correlation procedure. Templates may be constructed 
from both actual stellar spectra and lists of laboratory 
wavelengths. Since the former depend on spectrometer 
resolution, and the latter are subject to revision as bet- 
ter laboratory data become available, all such templates 
should be fully documented, as should the software used 
for the cross correlation (e.g., exactly what is being cor- 
related: the residual flux, or the line absorption; what is 
the weighting of different spectral orders; exactly how are 
th e observed line shifts co n verted into velocity values?). 

lOullberg fc Lindegren! l|2002l) deduced barycentric 
radial-velocity measures for some forty stars (using only 
Fei lines) with a median internal error of 27 m s _1 , and 
an external error of 120 m s _1 (the latter mainly com- 
ing from uncertainties of the wavelength scale in the solar 
spectral atlas used as wavelength reference). Although the 
precision achieved is somewhat lower than otherwise pos- 
sible, the accuracy is higher since the procedures involved 
are fully documented. Such radial- velocity measures there- 
fore become reproducible by other observers using dif- 
ferent instruments and different techniques, as evidenced 
by the good agreem ent for those stars in common with 
iNideveret all |2002l) . 

With improved measuring precisions, an increasing 
number of publications have started to use expressions 
of 'absolute' velocities, often meaning merely the use of a 
zero-point on the radial-velocity scale, obtained through 
calibrations against the solar spectrum or otherwise. The 
use of such a term is somewhat unfortunate since the con- 
cept of absolute velocity has a special physical meaning in 
relativity, denoting something rather more fundamental 
than, e.g., certain modes of calibrating wavelength-shift 
measurements. 



6.2.2. Data reduction and software 

The velocity values obtained as a result from spectroscopic 
observations depend not only on instrumental hardware 
effects, but increasingly also on the software versions used 
for reducing the data. 

Detailed procedures for the reduction of spectroscopic 
observations to the solar- s ystem barycentre have b e en de- 
velo ped e.g. by IStumpfj l|l977l Il979l Il98fl Il98.4 Il98fih 
and iMcCarthvl <) 19951) . based on th e series of solar-s ystem 
ephemerides available from JPL {Stan dishllT9901. The 
ephemerides and related computation services are con- 
veniently available on-line through JPL's HORIZONS 
System. 8 

At the 1994 IAU General Assembly it was decided 
to systematically set up software tools in order to en- 
hance the interchangeability of observational data and 
theoretical ideas. This set of tools is calle d the IAU 
Standards of Fundam ental Astronomy, SOFA I Fukushiinal 
1199.4 IWallacelll998h . The SOFA collection of algorithms 
include such for the accurate relativistic transformation of 



http:/ /ssd.jpl. nasa.gov/ 
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observed spectral-line displacements to the solar-system 
barycentre. 

Such algorithms apply to any periodic signals from a 
distant source, not only the periodic modulation inher- 
ent to an electromagnetic wave. In particular, the peri- 
odic modulation of pulsar signals follows the same math- 
ematical physics. In pulsar timing observations, the issues 
of uniform calibration of observations from different sta- 
tions, and their referral to the solar-system ba rycentre, 
have be en the topic of detailed examinations (Hellings 
Il986albl) . These transformations are for instance included 
in software packages such as TEMPO, a program for the 
analysis of pulsar timing data, maintained by Princeton 
University and the Australia Telescope National Facility. 9 
Many of these issues (including those of defining refer- 
ence frames, timescales, etc.) are directly applicable also 
to electromagnetic waves. 

The widesprea d FITS format (Flexible Image 
Transport System; IWells et alJ Il98l|) used for the han- 
dling of astronomical data is undergoing various modifi- 
cations and extensions, including a more elaborate rep - 
resentation of spectral quantities llOreisen et a .lJ l2nn,Sh . 
Alternative representations of spectral coordinates include 
'radio-convention velocities' computed from frequency 
shifts, 'optical-convention velocities' computed from wave- 
length shifts, 'relativistic Doppler velocities', and others 
(cf. Sect. I2.2jl . While the differences among such concepts 
may be small for most ordinary applications, any work 
aiming at very high accuracy should carefully examine the 
exact definitions of the various data fields, to understand 
how they can be transformed to barycentric radial- velocity 
measures. One has to remember that the prime purpose of 
standards such as FITS is not the accurate physical inter- 
pretation of data, but rather their transportation between 
different computers and software environments. 

6.3. 'Astrometric radial velocity': A stringent definition 
for geometric measurements 

The second resolution simply specifies how 'distance' and 
'time' should be defined in order to provide the geomet- 
ric measurement of radial motion called astrometric radial 
velocity. Briefly, the resolution states that the appropriate 
coordinate system is the Barycentric Celestial Coordinate 
System (BCRS, Sect. 14.10 . with time expressed as the 
barycentric time of light arrival f Sect. 1-0)1 on the barycen- 
tric coordinate time scale (TCB). Analogously, the con- 
ventional understanding of proper motion is generally un- 
derstood to mean the rate of change in barycentric direc- 
tion with respect to the barycentric time of light arrival, 
although we are not aware of any previous formal defini- 
tion to that effect. 



http:/ /pulsar. princeton.edu/tempo/ 



7. Unsolved issues 

The IAU resolutions were elaborated with the aim to 
permit results of spectroscopic and astrometric radial- 
velocity measurements to be unambiguously quantified on 
the 1ms -1 level. Many known effects on the sub-m s _1 
level are also taken care of in the present definitions. For 
example, the radial velocity of any object varies cyclically 
throughout the year, as the observer orbits the Sun and 
views the stellar velocity vector under a slightly differ- 
ent projection angle. Since the radial-velocity measure is 
defined relative to the solar-system barycentre, such ambi- 
guities are removed. However, there do exist other issues, 
where the present concepts may be inadequate. A few of 
them are highlighted below. 

7.1. Effects beyond the (inner) solar system 

The definition leaves 'uncorrected' all the (largely un- 
known) effects originating from outside the (inner) solar 
system. In particular, the BCRS describes an asymptot- 
ically flat metric at large distances from the Sun, thus 
ignoring effects of the gravitational fields from other indi- 
vidual stars and, on a larger scale, from the Milky Way 
Galaxy and structures therein (e.g., spiral arms, dark- 
matter concentrations). For instance, the large-scale grav- 
itational potential of the Galaxy causes wavelength shifts 
that may be relevant for highly accurate kinematic mod- 
elling. Over a range of galactocentric distances (i?) the 
galactic potential is cru dely described by that of a sin- 
gular isothermal sphere iBinnev fc Tremainelll987|) . lead- 
ing to a differential gravitational redshift between the 
star and observer of Av glav = (V 2 /c) ln(i? b s /-R*), where 
V ~ 220 km s _1 is the circular galactocentric speed. Thus, 
the spectra of stars in the central bulge (i?* ~ 1 kpc) 
may be gravitationally redshifted by 300-400 m s _1 , while 
stars in the Magellanic Clouds (i?» ~ 55 kpc) might be 
blueshifted by a similar amount, as seen by an observer 
near the solar position at i? b s — 8.5 kpc. 

7.2. Gravitational lensing 

Gravitational lensing, i.e., the bending or focusing of light 
during its propagation through gravitational fields may af- 
fect the radial velocity in different ways. In gravitational 
microlensing, another (fainter) star or other object passes 
very nearly in front of the target star (as seen by the ob- 
server) , and its gravitational field focuses light toward the 
observer. A stellar gravitational field is too weak to cause 
resolvable multiple images, so instead a source brighten- 
ing is observed. The velocities of stars in the Milky Way 
(acting as lenses) imply typical timescales for such events 
on the order of a few weeks (~ 10 6 seconds). 

The gravitational field of the lens causes a time delay 
of the light signal from the target star. This (Shapiro) de- 
lay, given by Eq. (|14() . is of order (20 /is) \n(r/p) for solar- 
mass lenses, where r is the distance to the lens (assumed 
to be half-way to the target) and p the impact param- 
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ctcr of the light ray. Thus for a ray grazing the stellar 
limb and observed at 1 kpc distance the delay is of order 
0.5 ms. A variation of the delay by this amount over a 
timescale of 10 6 s would cause an apparent change in the 
radial- velocity measure of the target star by ~ 0.15 m s _1 . 
Lensing by more massive or compact objects could thus in 
principle produce measurable variations. For a discussion 
of the corresponding effec t on pulsar timing observations, 
see lHosokawa et all 1 1999V Under certain conditions addi- 
tional relativistic effects causing time delays might enter, 
such as the Lense-Thirring or Kerr delay caused by the 
spin of the gravitational source and the ensuing frame- 
dragging. 

However, quite different effects could cause much more 
significant changes in the observed wavelengths of stellar 
spectral lines during a microlensing event. The amount of 
light amplification from the target object depends on the 
exact geometry of the target, the lens, and the observer. 
On this microarcsecond level, the disk of the target star 
is an extended object, and different parts of its disk grad- 
ually undergo different amounts of flux magnification, as 
the lensing object passes by. Since stars often rotate at a 
significant rate, portions of the stellar disk that approach 
the observer (with a spectrum Doppler-shifted to the blue) 
may at some time be differently enhanced from the red- 
shifted portions near the opposite stellar limb (receding 
from the observer) , producing a var iable wavelength shift 
on a level o f up to several km s _1 i Maoz fc Oouldlll994l 
iGouldlh 997ft . 

Gravitational lensing by more massive objects, e.g., 
clusters of galaxies, often produces multiple or extended 
images of the same target object. Each (sub)image cor- 
responds to a different light-path to the source, and thus 
a different Shapiro delay. For a geometry changing with 
time, there would be a variable differential delay, caus- 
ing each (sub)image to have a different (spectroscopic) 
radial velocity. Thus a particular source would not have 
one unique radial velocity, but different values depending 
on which among several light-paths from the source to the 
observer that are chosen. 



7.3. Gravitational waves 

The source of the gravitational lensing need not be sta- 
tionary, but could be transient, in the form of a passing 
gravitational wave. Although it appears that the effects 
will be very small, there might exist specific conditions 
(such as compact objects in close binary systems), where 
the variable time delays introduced by such waves s hould 
be tak e n into consideration, a t least in principle; e.g. , Fakir 
(|l994 : lKopeikin fe Ozernovl l|l999j) . 



motion relative to the general expansion, as represented 
by the Hubble parameter Hq __ 70 km s _1 Mpc -1 = 
70 m s _1 kpc -1 , or relative to the local expansion rate? 
Our accuracy aim of 1 m s -1 corresponds to the formal 
expansion velocity at a distance of only 14 pc, that of very 
nearby stars. 

The extent to which local systems participate in the 
general expansion of the Universe is a problem that has 
been t reated by several authors, beginning with lMcVittiel 
It has been argued that local entities such as the 
solar system or even the Milky Way Galaxy should be 
unaffected by the cosmic expansion since if everything ex- 
panded equally, the expansion would be unobservable. The 
full problem is quite complex, but there seems to be no 
fundamental reason why there should be a specific scale 
below which t here is no expansio n . For detailed discus- 
sions, see, e.g.. ICooperstock et al.l i|l998j) . and references 
therein. 

Finally, on cosmological scales of time and space, we 
cannot even be certain about the constancy of physical 
'constants'; to include such possible effects in the defini- 
tions remains a task for the future. 
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Appendix: IAU resolutions 

To enable high-accuracy studies of radial velocities, and 
to permit accurate comparisons between observers us- 
ing different methods, two resolutions were adopted by a 
number of Divisions and Committees of the International 
Astronomical Union, at a special session during its 
XXIVth General Assembly in Manchester (August 2000). 
The resolutions define a spectroscopic barycentric radial- 
velocity measure, and an astromet ric radial veloci ty. The 
full text of the resolutions follows l|Rickmanll200l . 
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Resolution CI on the Definition of a Spectroscopic 
"Barycentric Radial-Velocity Measure" 

Divisions I, IV, V, VI, VII, IX and X, and Commissions 
8, 27, 29, 30, 31, 33, 34 and 40 of the International 
Astronomical Union 



7.4. Cosmological effects 

A conceptual problem concerns the cosmological redshift: 
what is the meaning of 'radial velocity' in the context 
of an expanding Universe? Is it to be understood as a 



10 The texts of Resolutions CI and C2 are also available at 
http://www.astro.lu.se/~dainis/HTML/RADVEL.html and 
http://www.astro.lu.se/~dainis/HTML/ASTRVEL.html 
respectively. The resolutions referred to in the Note at the 
end of Resolution C2 are found in IAU Transactions XXIV 

B, pp. 37-43 and 44-49, and at 

http://danof.obspm.fr/IAU_resolutions/Resol-UAI.htm 
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Recognising 

1. That recently improved techniques for determining ra- 
dial velocities in stars and other objects, reaching and 
exceeding precision levels of meters per second, require 
the definition of "radial velocity" to be examined; 

2. That, due to relativistic effects, measurements be- 
ing made inside gravitational fields, and alternative 
choices of coordinate frames, the naive concept of ra- 
dial velocity being equal to the time derivative of dis- 
tance, becomes ambiguous at accuracy levels around 
100 m/s; 

Considering 

1. That, although many effects may influence the pre- 
cise shifts of spectroscopic wavelengths and frequen- 
cies, only local ones (i.e. arising within the solar sys- 
tem, and depending on the gravitational potential of 
the observer, and the observer's position and motion 
relative to the solar-system barycenter) can in general 
be reliably calculated; 

2. That, although the wavelength displacement (or fre- 
quency shift) corrected for such local effects can thus 
be derived from spectroscopic measurements, the re- 
sulting quantity cannot unambiguously be interpreted 
as a radial motion of the object; 

Therefore recommend 

That, whenever radial velocities are considered to a high 
accuracy, the spectroscopic result from a measurement of 
shifts in wavelength or frequency be given as the " barycen- 
tric radial- velocity measure" czb , after correcting for grav- 
itational effects caused by solar-system objects, and effects 
by the observer's displacement and motion relative to the 
solar-system barycenter. 

Here, c equals the conventional speed of light = 
299,792,458 m/s, and z B = (A - A )/A , where A is the 
rest-frame wavelength and A the wavelength observed by 
a hypothetical observer at zero gravitational potential, 
located at, and being at rest with respect to, the solar- 
system barycenter. The epoch of the observation equals 
the barycentric time of light arrival. 

The radial-velocity measure czb is expressed in velocity 
units: to first order in z-q it coincides with the classi- 
cal concept of "radial velocity", while avoiding the im- 
plicit interpretation as physical motion. The solar-system 
barycenter is defined by Resolution A4 adopted at the 
IAU XXIst General Assembly in 1991, and supplemented 
by Resolution B6 at the IAU XXIIIrd General Assembly 
in 1997. 



Resolution C2 on the Definition of "Astrometric 
Radial Velocity" 



Divisions I, IV, V, VI, VII, IX and X, and Commissions 
8, 27, 29, 30, 31, 33, 34 and 40 of the International 
Astronomical Union 

Recognising 

That recently improved astrometric techniques may per- 
mit the accurate determination of stellar radial velocities 
independent of spectroscopy, thus requiring a definition 
independent from spectroscopic measures; 

Considering 

That the change in the barycentric direction u to objects 
outside of the solar system is customarily expressed by 
the proper-motion vector /x = du/dte, where t B is the 
barycentric coordinate time (TCB) of light arrival at the 
solar system barycenter; 

Therefore recommend 

That the geometric concept of radial velocity be defined as 
dr/d^B, where r is the barycentric coordinate distance to 
the object and t B the barycentric coordinate time (TCB) 
for light arrival at the solar system barycenter. 

Note: The Barycentric Celestial Reference System (in- 
cluding the barycentric coordinate time) is defined in 
Resolutions B1.3 and B1.5 adopted at the IAU XXIVth 
General Assembly in 2000. 
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